Phases of anisotropic dipolar antiferromagnets 
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We study systems of classical magnetic dipoles on simple cubic lattices with dipolar and antifer- 
romagnetic exchange interactions. By analysis and Monte Carlo (MC) simulations, we find how the 
antiferromagnetic phases vary with uniaxial and fourfold anisotropy constants, C and D, as well 
as with exchange strength J. We pay special attention to the spin reorientation (SR) phase, and 
exhibit in detail the nature of its broken symmetries. By mean field theory and by MC, we also 
obtain the ratio of the higher ordering temperature to the SR transition temperature, and show 
that it depends mainly on D/C, and rather weakly on J. We find a reverse SR transition. 
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I. INTRODUCTION 

Long-range magnetic order brought about by purely 
dipolar interactions is somewhat rare in nature^ Inter- 
est in the subject is nevertheless growing. Some of it 
comes from the availability of synthesized crystals of 
organometallic molecules^ that behave at low tempera- 
ture as single spins Owing to the large organic mass 
enveloping the magnetic cores of these molecules, dipole- 
dipole interactions are then dominant £ Long range or- 
der has already been observed experimentallyi&SiiSiAiii^ 
Because uniaxial anisotropy is very large, Ising spins 
with dipolar interactions are reasonable models for these 
systems*^ Early rigorous work by Luttinger and Tisza 14 
established which type of magnetic order obtains at low 
temperature T in dipolar Ising models in each of the cu- 
bic lattices. The same results have been arrived at more 
recently by simpler methods^ 

No Ising model can however account for some 
of the interesting collective behavior, such as spin 
reorientationiSiiLiSiiSi 2 ^ (SR) and canted (i.e., non- 
collincar) spin configurations, that can be induced by 
purely dipolar interactions. Accordingly, three com- 
ponent spin models that include exchange as well 
as dipolar interactions are often used when modeling 
these effects To explore the transition between 

the dipolar- and exchange-dominated antiferromagnetic 
phases is one of the aims of this paper. 

The SR transition deserves special attention. Ther- 
mally driven SR transitions take place well within an or- 
dered magnetic phase. All spins rotate as a whole as the 
temperature decreases below some SR transition temper- 
ature T r . SR transitions can be continuous, as the ones 
first observed in the bulkiSiii or first order, as often ob- 
served in films liSii^iSSi It is important to realize that a SR 
phase comes with continuous transitions^ This phase is 
defined by its own unique set of broken symmetries. One 
of the aims of this paper is to exhibit this in detail. 

It has long been realized that higher order anisotropies, 
are required for the existence of the SR phase ,iL2i24 Ac- 



cordingly, we choose to study the model Hamiltonian 

H = Hj + H d + H A , (1) 

where TLj and H.^ are for all exchange and dipolar inter- 
actions, respectively, 

H A = -DJ2(SD 2 -C^iiSD 4 + (SI)% (2) 

i i 

D and C are constants, and Sf is the a component (a = 
x, y, z) of a unit spin S at site i. 
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FIG. 1: Phase diagram for T = 0. The type of order labeled 
c and s is for — l.SAcd < J < and J < — 1.34ed, respectively. 
Full and dashed thick lines stand for first- and second-order 
transitions, respectively. 

The fourfold anisotropy in Eq. (0 contributes the de- 
sired competition to the uniaxial anisotropy. Such a term 
arises naturally in tetragonal crystal lattices, and will be 
seen to bring about the canted spin phase that is shown 
in the lower left quadrant of Fig. 2] 

To continue with the definition of the model, we let, 
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where E(y) 1S a sum over an nearesi neighbor bonds, and 
Sj is a classical 3-component unit spin at lattice site i. 
We also consider nearest neighbor exchange interactions 
in order to establish how small | J | must be for dipolar 
magnetic order to obtain. 

Finally, H d = (1/2) Ey T^ P SfSf, where the 
sum is here over all sites i and j except i = j, 

* i 1 *-^) (4) 

and Tij is the displacement between sites i and j. We let 
Ea, Ej, and Ed be the values of Ha, H,j, and for 
whatever state is specified. 

Since we treat all spins as classical unit vectors, no 
quantum effect is taken into account. Accordingly, we 
disregard hyperfine interactions, since they have no ef- 
fect on the thermal equilibrium behavior of electronic 
spins when the latter are treated classically^ Quantum 
effects, such as the enhancement of the transverse field 
that is necessary in order to destroy long-range order in 
LiHoF4 below approximately 1 Kj2£ are therefore beyond 
the scope of this paper. Thus, some numbers we obtain, 
such as the crossover value of J, below which dipolar in- 
teractions become dominant, may change some when hy- 
perfine interactions are included in a quantum treatment. 
On the other hand, our main results, the inscnsitivity to 
the value of J to the ratio of the ordering temperature to 
the SR temperature, as well as the existence of a reverse 
SR transition, are of a rather qualitative nature, and are 
therefore expected to be insensitive to quantum effects. 

Simple cubic (sc) lattices and zero applied magnetic 
field H are assumed throughout. We only work with 
Lx Lx L box like systems, and let dipole-dipole interac- 
tions act between each spin and all other spins within an 
L x L x L box centered on it. We use periodic boundary 
conditions, because they give faster convergence towards 
the L — ► co limit than free boundary conditions do. In 
addition, because we work with antiferromagnets, and 
not ferromagnets, convergence is faster than would oth- 
erwise be. Most of our results follow from simulations for 
L = 8 and L = 16. This would be insufficient for critical 
behavior work but is adequate for our purposes, that is, 
establishing which type of magnetic order obtains and 
the corresponding phase boundaries. It is worth recall- 
ing that thermal equilibrium results obtained for H = 
for large cubic-shaped systems can, by virtue of Grif- 
fith's theorem^ be generalized to other shapes in three 
dimensions (3D). 

Our simulations follow the standard Metropolis Monte 
Carlo algorithm^ More specifically, we start simulations 
with an initial configuration in which all spins point in 
either random or parallel directions. We next compute 
the dipolar field at each site. Time evolution takes place 
as follows. A spin is chosen at random and temporar- 
ily pointed in a new random direction. Let AE be the 
corresponding energy change. If AE < 0, the tempo- 
rary direction becomes permanent. If on the other hand 



AE > the temporary direction becomes permanent 
only with probability exp(— AE/ksT), where T is the 
system's temperature. The field changes that ensue at 
every site in the system if the new spin direction is ac- 
cepted is then computed, thus updating the field values 
everywhere on the system. 

The plan of the paper and a list of the results ob- 
tained follow. In Sec. [HJ wc study all the phases, except 
the SR phase (see Fig. 0, for all J < 0. To this end 
we first define a canted state which lets all spins point 
along the easy magnetization axes for any C and D while 
Ed remains invariant if the system's shape is cubic. This 
is the main device which enables us to conclude in Sec. 
ITT1 that c and s, exhibited in Fig. 2] are shown to be 
the ground state configurations for — 1.34e<j < J < 
and J < — 1.34ed, respectively. We also report MC data 
which suggests that, as in the Landau theory^ the xy- 
collinear to a;y-canted transition line is first order. We 
also report results from MC simulations for temperature 
driven phase transitions between paramagnetic and xy- 
collinear and between xy-collinear and xy-canted phases. 
In Sec. IIIII we study the SR phase. In an anisotropy di- 
agram we show the free energy minima that follow from 
MC simulations for the SR phase. The broken symme- 
tries are clearly exhibited. This suggests how to write 
mean field equations for the SR transition, from which 
we obtain T r as well as the higher ordering temperature 
as a funtion of D/C. We show that the ratio of these 
two temperatures depends mainly on the D/C ratio and 
rather weakly on J, and mean field yields approximately 
the same value for it as MC simulations do. For D < 
and 0.8 < D/C < 1, cooling through the paramagnetic 
phase first brings the system into xy-canted ordering. 
The SR phase is encountered at a lower temperature. For 
D < but < D/C < 0.8, cooling through the para- 
magnetic phase first brings the system into z-collinear 
ordering, and the SR phase is encountered at a lower 
temperature. An interesting phase diagram obtains at 
D/C ~ 0.8, and a reverse spin reorientation that is tem- 
perature driven is observed. A complete spin reorienta- 
tion, from the easy magnetization axis into the perpen- 
dicular plane, first takes place upon cooling below the 
paramagnetic phase, followed, upon further cooling, by a 
reverse spin reorientation towards the easy magnetization 
axis. This is so for J = as well as for Ed = 0. 



II. COLLINEAR AND CANTED PHASES 

In this section we study all the phases shown in Fig. 
n except the SR phase. We first study how the ground 
state of a purely dipolar system varies with D and C . 
We then extend this to all J < 0, and finally consider 
nonzero temperatures. 



3 



A. The canted state 

In this subsection we define the canted state c. We do 
this in order to avoid the difficulty which follows from 
the the fact that dipolar interactions are not rotationally 
invariant, thus making the task of simultaneously mini- 
mizing Ea and E& nontrivial. The canted state c enables 
one to do this minimization. Let 

Sf = rf cos 9, Sf = t? sin 9 sin cf>, Sf = 7? sin 9 cos 0, 

(5) 

where 9 is the angle between the spin vector and the z- 
axis, <j) is the azimuthal angle, Tj = [rf, rf , r*] is given 
by 

T . = u_\\y(i)+z{i)^ (—i^M+vtfh (6) 

[x(i), y(i), z(i)] is the position of site i and x(i), y(i), and 
z(i) are all integers. 

Before proceeding with the argument, we point out 
some basic features of the canted state we have just de- 
fined. Note how Ti varies with i. A canted spin config- 
uration therefore follows from Eqs. (J5J and ijfjj if 9 ^ 
and 7^ 7r/2. For — ir/2, a canted spin configuration 
also obtains if 4> ^ and </> ^ 7r/2, as, for instance, the 
c state in the lower left hand quadrant of Fig. 2] The 
fact that our definition of a c state gives a collinear spin 
configuration for these few special cases (see the c state 
in the lower right hand corner of Fig. Q]for 9 = 0) should 
not be too confusing. 

To start our argument, note that 9 = gives the 
minimum value of E& for sc lattices in the bulkiii2£ 
It is next shown that the dipolar energy E& of cubic- 
shaped systems is independent of 9 and of <fi in the 
canted stated In order to see this, consider first Sj • Sj, 
which upon substitution of Eq. (J5J becomes t/tJ cos 2 9+ 
rfrj sin 9 sin <f> + rfr x sin 2 9 cos 2 <j>. Now, it is easy to 
check, that the sum ^« ■ Sj/(r,*j), where f(tij) is 
any function of ry, over a sc lattice bounded by the 
surface of a cube is independent of 9 and of 0, since 
E^t/t//^-) = E^rM/fa) = E^rfr;/(ry) 
then. Thus, the first of the two terms of the dipo- 
lar interaction is independent of 9 and of cf>, and, by 
the way, so is the exchange energy. Now, consider 
Wij = (xijSf + HijSf + ZijSf){xijS^ + VijSj + ZijS*), 
from the second term of the dipolar interaction. Again, 
w ij/ r ij over a sc lattice bounded by the surface of a 
cube is independent of 9 and of <f>, since (1) cross terms 
do not contribute to the sum, by reflection symmetry, 
^d (2) T i T j/ r ij = = EyTfrf/r|. 

Therefore, 9 and <f> can be freely chosen in Eq. (J5J in 
order to minimize Ea, at no cost to Ed, which is the 
desired result for bulk systems. 

Finally, since the minimum value of Ea can be reached 
by appropriate choice of 9 and of </> in Eq. JSJ) , it follows 
that minimization of Ea with the c state gives the ground 
state. 

The phases obtained for the ground state by this pro- 
cedure if J = are shown in Fig. ^ We refer to the 
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FIG. 2: (Color online) Order parameter m™, for a = x,y,z, 
vs C for systems of L x L x L dipoles on sc lattices at T = 
0.04£d/fcs, where fcs is Boltzmann's constant, for D = — 3e<i. 
A transition between the a;i/-collinear and xy-canted phases 
is clearly exhibited at C = 0. In order to reach equilibrium, 
T was lowered from T = 2ed/kB in AT = — 0.02ed/kB steps. 
2 x 10 4 MC sweeps were made at each temperature step. Equi- 
librium values come from averages over 2 x 10 5 MC sweeps. 



states on the right hand side, upper left hand side and 
lower left hand sides of the diagram as z-collinear, xy- 
collinear, and xy-canted phases. The nature of the SR 
phase is the subject of Sec. IIIII 

Note that the xy-collinear phase is unstable with re- 
spect to any nonvanishingly small C < 0. [On the other 
hand (see below), the xy-canted phase obtains in finite 
nonzero temperature only for a sufficiently small C < 0.] 

Before we proceed any further, we need some addi- 
tional definitions. Let 

m^N-^Sfr?. (7) 

i 

In a c state, m\ = cos6>, m x t = sin 9 cos cf>, and rwj. = 
sin 9 sin <fi. 

According to Landau's theory, the transition between 
the z-collinear and the xy-plane phases, as well as the 
transition between the rry-collinear and the xy-canted 
phases, is of first order m& This has been observed for the 
z-collinear to xy-plane transition in films. 19 - 22 30 The re- 
sults we have obtained from MC simulations also support 
this conclusion. 

On the other hand, we know of no experimental or 
MC evidence about the nature of the xy-collinear to xy- 
canted transition. Data points from MC simulations we 
have performed are plotted in Fig. [21 Note: (1) = 
for all -1 < C/e d < 1, (2) m x z = m y c for all C/e d < 0, 
and (3) either (a) m x = and m v t = 1 or (b) m x = 1 
and m v c = for C > 0. A first order phase transition, 
between two in plane spin configurations, one along the 
x or y axes and the other one with <fi = tt/2, is clearly 
suggeted. 



B. Exchange interaction 

The effect exchange interactions have on the ground 
state is studied in this subsection. 

Let's first define the collinear s state, illustrated in Fig. 
Hby one of the two sets of the states. In this state, all 
nearest neighbor spins to spin S point opposite to S (as 
in Fig. Alternatively, in a s state, 

Sf = rji cos 6, Sf = rji sin 9 sin cj), Sf — f]i sin 9 cos (f>, 

(8) 

where 

T)i = (— l^W+fW+^O (9) 
It makes sense to also define 

m? = iV- 1 ^Sfr ?i . (10) 

i 

In a s state, = cos 9, nig = sin 9 cos </>, and mi = 
sin 9 sin 0, but in a c state, m s = 0. 

We showed in Sec. Ill Al that c is the ground state if 
J = and that Ed is independent of 9 and <fi in Eq. 
(0. In addition, E d = -2.67e d theni^ Note also that 
Eqs. JSJ and © imply Ej = J for all 9 and in a 
c state. Similarly, it can be easily checked that in a 5 
state, Ed — 0, independently of the direction of S, and 
that, Ej = -3J. 

Therefore, of the above two states, the c state gives 
the lower energy in the — 1.34e ( ; < J < range and the 
5 state gives the lower energy in the J < — 1.34£d range. 
Monte Carlo simulations show that no other state gives 
a lower energy in the whole J < range, that is, 

E = -2.67s d + J + E A (11) 

if -1.34e d < J < and 

E = 3J + E A (12) 

if J < -1.34g d . 

Clearly dE/dJ is discontinuous at J ~ — 1.34e<i, im- 
plying, a hrst order transition betwwen purely dipolar 
induced canted (i.e., a c state) and collinear antiferro- 
magnetic phases (i.e., a s state) for all C and D. 

C. T / 

Of all thermally driven transitions between any two 
of the paramagnetic, xy-collinear, and xy-canted phases 
(the SR phase is treated in Sec. IHIfl . the most interesting 
one is the one between the xy-collinear and the ccy-canted 
phases. It is illustrated in Figs. and 0] This effect 
has been discussed before for thin films^i Note how the 
variation of | [mf] 2 — [m^] 2 | and of | m x c m v c \ at T ~ 
Q.2Aed becomes sharper as the systems size increses. This 
suggests a first order phase transition, in accordance with 
Landau's theory^ 
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-1.5 -1.0 -0.5 0.0 0.5 1.0 

-C/D 

FIG. 3: (Color online) Transition temperatures and phases 
for systems with only dipolar interactions (J = 0), and two 
values of D: —3Ed (A and 0), and —0.5^ (o and □). + are for 
the SR phase boundary. The data come from MC simulations 
of systems of 8 x 8 x 8 spins on sc lattices. The vertical dashed 
line stands for the cooling path that was taken to obtain the 
data points shown in Fig. 




T [e d /k B ] 

FIG. 4: (Color online) | \rn c f ~[m v c } 2 \ [2m^m v c ] 2 and C/k B vs 
T for systems of L x L x L spins with J = 0, and D — — 0.5e<;, 
and C — 0.05D. For L — 8, + (red online), x (red online), 
and A (red online) stand for | [m*] 2 - [m y c ] 2 |, [2m^m?f 
and C/kB, respectively; for L = 16, □, 0, and o stand for | 
[m x c \ 2 -[m y c \ 2 |, [2m>?] 2 and C/k B , respectively. For L = 16 
(L — 8), all data points stand for averages over 10 5 (2.5 x 10 5 ) 
MC sweeps. 



III. THE SPIN REORIENTATION PHASE 

A. T = 

We first explore SR at T = 0. We start by minimizing 
the anisotropy energy Ea for all C and D in the SR 
phase (see Fig. Recall from Sec. [H]that Ea can be 
minimized freely for all J < 0. We obtain 

S = (±u,±u,±v) (13) 
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-D/C 

FIG. 5: mf vs —D/C for a = x, y, z and T = 0. • stands for 
both mf and rrv/, and ■ stands for to*, for C = — 0.5ed and 
J = 0; □ stands for both mf and m\ , and o stands for mf , for 
£d = and C = J < 0. All symbols are from MC simulations 
in which T was lowered from T 3> £<j, | J | to T -C e<j, | J |. 
Continuous lines are from Eq. 11311 and below. 

for C < L> < 0, where u = y/D/2C and v = y/l - 2u 2 . 
Thus, varying D through the C < D < range forces 
the spin directions to vary in the ground state from the 
z-collinear state shown in Fig. to the xy-phase. This 
is analogous to the phenomenological theory of SR^ 

We have simulated cooling from high temperatures 
down to ksT <C Ed, J for various values of C and D: (1) 
purely dipolar interacting 3D systems (i.e., J = 0); (2) 
antifcrromagnetic 3D systems with only nearest neigh- 
bor interactions (i.e., J < 0, Ed = 0). The MC results 
obtained as well as the numbers that follow from Eq. 
(|13|l are plotted in Fig. \5\ The data points fall on the 
predicted curves, independently of J. 

Note that Eq. (|13|l allows 8 different spin directions. 
This suggests which symmetries are broken in the SR 
phase. We examine this idea in more detail in the follow- 
ing section. 

B. The temperature driven transition 

In Fig. mf for a = x,y, z, and the specific heat C 
are plotted vs T for J = 0, and D = — 0.5e f j, C = —Ed- 
The data points are from MC simulated cooling. Data 
points obtained from heating from ordered states, using 
Eq. (|13|) and below, do not differ from the ones shown. 

Specific heat curves obtained for various system sizes 
are shown in Fig. The data are not good enough for an 
accurate value of any critical index, but they suggest the 
transitions are continuous at both of the two transition 
temperatures. 

We have obtained additional specific heat and m c 
curves (as the ones in Fig. for J = and C = —3£d, 
from further MC simulations. The data points shown in 
Fig. [HI for the two transition temperatures as well as the 
labels shown for the phases follow from such curves. 




T [e Ik ] 

FIG. 6: ml (♦), m\ (■), to? (□), and C/k B (o) vs T. All 
data points come from MC simulations of 8 x 8 x 8 spins on 
sc lattices for J — 0, C = —Ed, and D = 0.5C. At each value 
of T, at least 10 s MC sweeps were made. Lines are guides to 
the eye. 




T [e d /k B ] 

FIG. 7: Specific heat vs temperature for systems of L x L x L 
dipoles, with J = 0, C — —Ed and D = 0.5(7. o, □ and o 
stand for L — 16, 8 and 4 respectively. T was lowered in AT = 
—0.01 steps. At each value of T, we made 5x 10 4 ,5x 10 s and 
2 x 10 7 MC sweeps for L — 16, 8 and 4, respectively. The high 
(low) temperature peak corresponds to the paramagnetic to 
z-collinear phase (z-collinear to SR phase) transition. Inset: 
log Cmax/ks vs log L for both transitions, o (□) stands for 
the high (low) T peaks. The slope of the dashed line is for a 
first order transition (C ~ N then). The slope of the full line 
is 0.2, in accordance with a countinuous transition in which 
a/v w 0.2. 



For all J < 0, the phase diagram shown in Fig. [SJis 
qualitatively the same, except that for J < —l.SEd, the 
xy— canted phase is replaced by the collinear xy phase 
shown in Fig. \T\for C,D < 0. 
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FIG. 8: (a) Phases of dipolar (J = 0) antiferromagnets for 
C = — 3ed- All data points come from MC simulations of 
8x8x8 spins. The dashed line is for the cooling path from 
which the data points in Fig. HJ follow. 



C. Broken symmetries 

In this section we illustrate the symmetries of the var- 
ious magnetic phases a system goes through when cool- 
ing from the paramagnetic into the SR phase. To this 
effect, a statistical sample of values m c takes up through- 
out time in each of the three phases -paramagnetic, 
z-collinear, and SR— that obtain for C = —Ed and 
D = 0.5C is shown in Figs. GJa) and^b). 

Now 



?K) = z 



(14) 



where p(m c ) is the probability to find a m c value 
and F(m t ) is the free energy. It follows that cloud 
densities in Figs. EJa) and Etb) are proportional to 
cxp[— i 7, (m c )/fcsT]. The clouds therefore stand for neigh- 
borhoods of the minima of F(m t ) for each of the three 
phases. 

The data points shown in Fig. EJa) for each value of 
T come from one single MC run. This is in contrast to 
the procedure used to obtain the data points shown in 
Fig. E[b), which were obtained from 10 3 independent 
MC runs of 4 x 10 4 MCS steps each. 

Clearly, all the symmetries that are broken in the xy- 
canted phase, in addition to the symmetries that are bro- 
ken in the z-collinear phase, are broken in the SR phase. 



D. Mean field 

It is interesting to calculate the mean field transition 
temperatures associated with the SR phase. We assume 
a magnetic order as either in a canted or a collinear state, 
and that only m c and m 5 remain to be determined. More 
precisely, we assume 



{ST) 



(15) 



(a) 




(b) 




mj 



FIG. 9: (Color online) (a) m c values measured every 10 
equally spaced MC sweeps in equilibrium for systems of 
8x8x8 dipoles with J = 0, C = —Ed and D = 0.5C. Each set 
of o, gray (red online), and black dots stand for ksT/ed = 0.9, 
0.5, and 0.15, respectively. All data points for each T come 
from a single MC run of 4 x 10 6 MC sweeps each. In each MC 
run, before we take any data, we let the system equilibrate 
by lowering the temperature from T = 2ed/kB (deep in the 
paramagnetic phase) in steps of AT = —0.1, of 4 x 10 3 MC 
sweeps each, (b) Same as in (a) but each of the three sets 
of points comes from 10 3 independent short MC runs (i.e., 
first cooling from T = 2ed/kB in each of the 10 3 MC runs) of 
4 x 10 4 MCS steps each. 



We now write the mean field equations. First note that 
multiplying the above equation by r" and summing over 
two nearest neighbor sites along the a direction gives, 



2m? = ]T t?{S? 



Similarly, one obtains 



2<= E ViiS* 

Q,i=l,2 



(16) 



(17) 
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where 



FIG. 10: Mean field transition temperatures T z and T xy and 
the nature of the antiferromagnetic phases vs —D/C for 3D, 
nearest neighbor interactions only (i.e., Ed — 0), J < 0, and 
C = 47. 



We now write the main mean field equation, 

(Si) = Z-i^Tuexph^uV^T] (18) 

u 

where the sum is over all directions of the unit vector u, 
e i (u) = -h(i)-u + e i4 (u), (19) 



e A (u) = -Du 2 z -C(ui + u 4 v ), 
h(i) = hj(i) +h d (i), 

h a J (i) = J{m«Ys T i+<Y. r )il 

ni ni 

the sum here is over all nearest neighbors of site z, 



j 



and, finally, 



Z = ^exp[- £l (u)/fc B T]. 



(20) 
(21) 

(22) 
(23) 



z= dn e - £A(u K 



(25) 



and J d£l is an integral over all directions of u. Simi- 
lar equations obtain for the the paramagnetic-xy-canted 
phase boundary. For other phase boundaries, lineariza- 
tion cannot be carried out for all spin components, 
whence slightly more complicated systems of equations 
ensue. We have solved the equations for the phase 
boundaries numerically. The solutions for Sd — and 
C = 4 J < are shown in Fig. EH 
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FIG. 11: (Color online) Temperature ratio T xy /T z vs -D/C, 
for the shown values of Ed, J, and C in cubic-shape systems. 
MF stands for mean field theory results. • and ▲ stand for 
8x8x8 and 16 x 16 x 16 spin systems, respectively. Data 
points from MC simulations are for 8x8x8 spins. Lines are 
guides to the eye. 

It is interesting to compare the ratio T z /T xy between 
the two relevant temperatures for the SR transition. 
Some results we have obtained from MC simulations and 
from the mean field equations above are shown in Fig. 
im How little the results seem to depend on whether SR 
is brought about purely by dipolar interactions (J = 0) 
or solely by exchange interactions (ea — 0) is noteworthy. 
This may be related to the fact that SR is completely in- 
dependent of J in the ground state (see Sec. Hljl . How 
close the mean field and MC data points lie over most of 
the < D/C < 1 range in Fig. E3is a l so intriguing. 



The above equations simplify near the phase bound- 
aries. Consider, for instance, the boundary between the 
paramagnetic and the z-collinear phases in Fig. 1101 Let 
T z define the the boundary between a phase where ml 
and m z s vanish and a phase where either ra\ or does 
not. T xy is similarly defined for the x and y compnents 
of m 5 and m c . 

Then, exp(— h • u^/feeT) can be linearized, and 



(24) 



E. Reverse SR 

The rather abrupt variation of T z /T xy near D/C = 0.8 
(see Figs. l8*land HT|l is worth exploring. 

Consider first J = 0. The behavior that obtains as 
a function of T for D = — 0.747C, and, somewhat ir- 
relevantly, C = —3ed, is illustrated in Figs. IT2T a') and 
Efb). Note that (| m\ |) seems to approach zero in the 
0.35 < T < 0.55 Ed/ks range as system size increases. 
Thus, four phases are visited as T varies in Fig. I12t ah 
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FIG. 12: (a) mf , for a = y,z vs T for systems of L x L x L 
dipoles with J = 0, C — —3ed and D = 0.747C. Data points 
for ml would fall right on top of data points for and are 
not shown. Open (closed) symbols stand for data points from 
MC simulations in which T was lowered (raised), starting from 
T — Ed (T = O.Oled). Circles and rhombi stand for a = z and 
x, respectively, both for L = 16. + (x) stand for a = z and 
a — y for systems with L — 8. All data points come from 
averages of | m" | over some 10 5 MC sweeps, (b) Same as in 
(a) but for C/k B vs T. 



the paramagnetic phase, at 0.7ed/kB ^ T, the z-collincar 
phase, at 0.55 < T < 0.7 Ed/kB, the rcy-canted phase, at 
0-35 < T < 0.55 e d /k B , and the SR phase, at T < 0.35 
Sd/ks- We refer to the SR that takes place within the SR 
phase, where T < 0.35e<j/^_B in Fig. I12f a). as a reverse 
SR, since spin orientations then vary opposite to the way 
they vary at the higher transition (near T = 0.55ed/kB)- 

The rather abrupt change in ttv( and in m z c near 
T = 0.7 is in accordance with a first order transition, as in 
Landau's theory. 23 (This is so in Landau's theory because 
neither symmetry group in either of the two phases is a 
subgroup of the other one. Note that the opposite con- 
dition obtains at the other two transition points, hence 
continuous transitions are predicted therein.) 

Reverse SR does not seem to be an isolated phe- 
nomenon. It also obtains for other values of D, in the 
range 0.72 < D/C < 0.77, and for nonzero values of the 
exchange constant. This is illustrated in Figs. I13f a) and 
H3T bL for Ed = 0, D/C = 0.8, and, less importantly, 
C = 4J. From similar plots to the one shown in Figs. 
I12f a) and !12f b). we have obtained for various values of 
C/D the phase diagram shown in Fig. ITU 
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FIG. 13: (Color online) (a) m" for a = y, z vs temperature 
for systems of L x L x L spins with Sd — 0, C = 4J < 0, and 
D = 0.8C. Data points for fall on top of the data points 
for rrij and are therrefore not shown, m" is represented by 
(red online) o, +, and continuous line for L = 32, 16, and 8, 
respectively, mj is represented by (black online) (}, x, and 
dashed line for L — 32, 16, and 8, respectively. All data points 
come from averages of | mf | over 1.8 x 10 J MC sweeps, (b) 
Same as in (a) but for C/kB vs T. 
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FIG. 14: Magnetic phase boundaries for systems of L x L x L 
spins in which J — and C = — 3ed- The full (dashed) line 
and + (o) is for L = 16 (L = 8). The vertical dashed line with 
an arrow is for the cooling path from which the data points 
in Fig. [TU follow. 



IV. CONCLUSIONS 

We have studied the effect of a fourfold anisotropy 
on the magnetic phases of systems of classical mag- 
netic dipoles on simple cubic lattices with dipolar and 
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antiferromagnetic exchange interactions. For negative 
anisotropy constants D and C , and for T = 0, we find 
canted and collinear spin configurations, as shown in Fig. 
[T] if -1.34ed < J < and J < -1.34s d , respectively. An 
interesting temperature driven transition between two or- 
dered phases is also reported (see Fig. 0J. 

We have studied in some depth the spin reorientation 
phase. Its broken symmetries are exhibited in detail. 
In our model, a thermally driven continuous SR transi- 
tion occurs if, in the ground state, spins are tilted away 
from the crystalline axes. We show, by mean field the- 
ory and by MC simulations, that, upon cooling below 
the paramagnetic phase, collinear (along the easy mag- 
netization axis) spin configurations obtain if D < and 
< D/C < 0.8, and spin reorientation (towards the per- 
pendicular plane) takes place at lower temperatures. On 
the other hand, if if D < and 0.8 < D/C < 1, spins 
point perpendicularly to the easy magnetization axis be- 
low the paramagnetic phase, and a spin reorientation to- 
wards the easy magnetization axis takes place below a 
lower transition temperature. 

Intriguingly, the ratio between the two relevant tran- 



sition temperatures, T z /T xy (see Fig. ITT)l . seems to de- 
pend, over most of the < D/C < 1 range, neither on 
the strength of the exchange interaction, nor on C, nor 
on whether the two temperatures follow from MC simu- 
lations or from mean field thoery. 

If D < and D/C ~ 0.8, then (see Figs. [13 U3 and 
I14fl . independently of the strength of the exchange inter- 
action, a complete spin reorientation, from the easy mag- 
netization axis into the perpendicular plane first takes 
place upon cooling below the paramagnetic phase, fol- 
lowed at lower temperature by a reverse spin reorienta- 
tion, in the SR phase, towards the easy magnetization 
axis. 
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